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Abstract: Experiments searching for weak interacting massive particles with noble gases such 
as liquid argon require very low detection thresholds for nuclear recoils. A determination of the 
scintillation efficiency is crucial to quantify the response of the detector at low energy. We report 
the results obtained with a small liquid argon cell using a monoenergetic neutron beam produced 
by a deuterium-deuterium fusion source. The light yield relative to electrons was measured for six 
argon recoil energies between 11 and 120 keV at zero electric drift field. 
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1. Introduction 

It is known from astronomical observations that most of the gravitational mass in the universe 
is made of dark energy and non-baryonic dark matter which does not couple to electromagnetic 
radiation [HI]- Dark matter has survived since the birth of the universe and hence must be stable 
and weakly interacting. The most prominent candidates for dark matter are the Weakly Interacting 
Massive Particles (WIMPs) [0], in particular the spin-1/2 neutralino with predicted mass in the 
GeV to TeV range. The neutralino would scatter on constituent quarks in nucleons, leading to 
nuclear recoils in the range of 1 - 100 keV. Non-accelerator laboratory (“direct”) searches are all 
based on the detection of nuclear recoils. The scattering cross section on nucleons is tiny, in the 
range 10“^ to 10“^^ pb, comparable to that for neutrino interaction. For low masses the sensitivity 
decreases due to the low recoil energy and the detection threshold, while for high masses the loss 
of sensitivity is due to the exponentially diminishing WIMP flux. 

The laboratory observation of dark matter is one of the most pressing issues in Particle Physics. 
The neutralino is being searched for at the LHC and direct searches for WIMPs are underway in 
non-accelerator underground experiments (for a review see [0]). The most stringent upper limits 
for WIMPs are obtained for the spin independent WIMP-nucleon cross section. Upper limits on the 
WIMP-nucleon cross section have been obtained by the XENONIOO 0] and LUX [0] experiments 
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using liquid xenon, the latter quoting the most stringent upper limit of 7.6 x 10“^^ pb for a WIMP 
mass of about 33 GeV. 

At the LHC the neutralino would manifest itself by leaving a large missing energy when scat¬ 
tering off quarks. These “indirect” searches at the LHC are more competitive than direct searches 
for low WIMP masses, providing the best upper limits of about 10“^ pb below 3.5 GeV [@, 0]. 

Several direct searches were performed or initiated during the last ten years, employing dif¬ 
ferent types of detectors such as solid state detectors or noble liquids (argon, neon or xenon) [ 0 ] 
which can be used in large volume time projection chambers. Experiments using liquid argon are 
underway: DEAP [H] and MiniCLEAN [0 at SNOLAB, DarkSide [110] at Gran Sasso, and ArDM 
firrH at Canfranc. 


2. Liquid argon for direct searches 


In liquid argon (LAr) scintillation light is emitted in a narrow VUV band around 128 nm with two 
components of different decay times, Ti ~ 7 ns and T 2 ^ 1.6 /is [1T211 . Heavily ionizing projectiles 
such as a particles or nuclear recoils contribute mostly to the fast decaying component, while the 
contribution from electrons and 7 -rays to the slow component is larger. In this work we define 
the prompt fraction Fp as the fraction of prompt over total integrated pulse height, the prompt 
light being chosen as the integrated pulse height in the time interval between —20 and -h30 ns of 
the time at which the light pulse reaches its maximum. Eigure shows the prompt fraction Fp 
measured during the early stage of the experiment, when a ^^^Po a-source (5.3 MeV) was installed 
is typically 0.25 and 0.75 for electron recoils and nuclear recoils, respectively. The prompt 
fraction is therefore useful as a discriminant to identify nuclear recoils and to reject background 
flU. This property is the main advantage of argon over xenon (in addition to its lower price), but 
the disadvantage is the presence of the radioactive j3-emitter ^^Ar. 
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Figure 1. Fraction Fp of prompt to total light for various projectiles. 


The scintillation light yield and its dependence on nuclear recoil energy is an essential input 
to the WIMP-nucleon cross section. The relative scintillation efficiency Le/f is defined as the ratio 
of the scintillation yield L„r from nuclear recoils to that from electron recoils L^r, 
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where T^r and Ter stand for the nuclear recoil and electron recoil energy, respectively. Ler is deter¬ 
mined from photopeaks or Compton edges, using y-sources. The relative efficiency Leff being an 
energy dependent quantity, a standard reference for calibration is required. For the sake of conve¬ 
nience we use in this work the 60 keV electronic photopeak from a source. As discussed 

below, we have measured the light yield from various y-sources (including a 32 keV ^^^Kr source) 
and have established its linearity for recoil electrons above 30 keV. The results obtained here can 
therefore be compared directly with those obtained from earlier experiments [O |T3,0] using e.g. 
the Ter = 122 keV photopeak from a ^^Co source OH]. 

The relative scintillation efficiency is measured by mimicking the WIMP-nucleus interaction 
with neutron beams, e.g. by elastic scattering of monoenergetic neutrons under fixed scattering 
angles (see QH, O] for argon and flTSII for xenon data). The measured light yield is then compared 
with a simulation and Leff is derived by applying iterative fits. This article reports the determination 
of Leff in LAr at zero electric field using a monoenergetic neutron beam from deuterium-deuterium 
fusion. Details can be found in flT3l1 . The measurements were motivated by the envisaged DARWIN 
dark matter project fEOl . 

2.1 Experimental setup 

A good way to calibrate the light output and to study the response of LAr to nuclear recoils is 
to scatter a beam of monoenergetic neutrons of energy T^ off argon nuclei, and to measure the 
light yield as a function of scattering angle 0, from which the recoil energy T^r can be calculated 
according to the formula 


( 2 . 2 ) 

where © is the scattered neutron angle in the laboratory and A 1 is the atomic mass number. The 
method is illustrated in figure ^ (left). 




Figure 2. Left: principle of the scattering experiment with monoenergetic neutrons scattered under the angle 
0. Right: nuclear recoil energies in argon for incident neutron energies of 2.45 MeV. 


We are using monoenergetic neutrons of T^ = 2.45 MeV from an electrostatic fusion source 
based on the reaction dd ^^He n. The energy deposits are plotted in figure § (right) as a function 
of scattering angle. The commercially available neutron generator is of the deuterium-deuterium 
plasma fusion type from NSD-fusion [ETII . The source is a cylinder at ground voltage containing 
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low pressure deuterium gas (10“^ mbar). An internal perforated cylindrical electrode at high volt¬ 
age (typically 80 kV) induces a discharge. The ionized deuterons are accelerated towards the inner 
electrode and accumulate in the central regions with energies of about 15 keV. This is sufficient to 
overcome the Coulomb barrier and induce fusion. The plasma of ^25 mm length emits neutrons 
isotropically into An with a typical rate of 10^ n/s. The certified operation time is 25 000 hours, 
after which the cell must be exchanged. The fusion generator is surrounded by a 90 cm diameter 
shield of borated polyester and the experiment is confined within a radiation controlled fence in 
our laboratory at CERN (figure ^ left). Residual radiation (mainly from scattered neutrons and 
X-rays) is well below the authorized limit of 2.5 /iSv/h. 



Figure 3. Left: neutron gun, argon cell and liquid scintillation counter used to measure Leff in LAr. Right: 
sketch of the LAr cell with its vacuum chamber (see text). 


The neutrons are collimated through a polyethylene orifice within roughly 0.2% x An sr to 
match the sensitive volume of the LAr cell. The collimator is defined by a lead box (wall thickness 
2 mm) inserted into the polyester to protect against X-rays. The LAr cell (figure right) is located 
at a distance of 73 cm from the neutron emission point. The sensitive volume of the cell consists 
of a 1.5 mm thin aluminium cylinder (75 mm in diameter and 47 mm high, hence a LAr volume of 
0.2^). The cell itself is contained in a larger cylindrical vessel (141 mm in diameter), also filled with 
LAr. The inner wall of the cell is covered by a Tetratex foil coated with a wavelength shifter made 
of tetraphenylbutadien (TPB) to shift the VUV scintillation light to a longer wavelength, following 
our developments described in fE2[1 . The surface density of the TPB covering the reflector is 1.0 
zb 0.1 mg/cm^. The light readout is performed with two Hamamatsu R6091-01 photomultipliers 
(PMT) on each side of the cell. The 3” bialkali PMTs have 12 dynodes and a Pt underlay to reduce 
the resistivity of the photocathode when operating at LAr temperature (^86 K). The surfaces of the 
PMTs are coated with a TPB/paraloid mixture (quantum efficiency estimated to be around 15%). 
To perform measurement at zero electric field the aluminium cylinder is polarized at the same 
voltage as the photocathode of the PMTs. 
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The detector is filled with condensed pure argon gas 6.0 containing less than 0.5 ppm of H 2 O, 
0.1 ppm of O 2 , 0.1 ppm of H 2 , 0.1 ppm of CO 2 , and 0.1 ppm of CO. Platinum resistors are used 
to monitor the LAr level inside the cell. The vacuum and gas pressure are read out by pressure 
sensors. A blue LED (A = 390 nm) connected to an optical fiber is inserted into the vessel for 
gain calibration. A slow control program written in Labview is used to monitor and record the 
temperature and pressure during operation. 

For continuous operation a gas system is connected to the cell to provide condensation and 
recirculation of the argon. The gas is condensed on the top of the chamber by the cold head of a 
Sumitomo CH210 cryocooler with a nominal cooling power of about 80 W, driven by a Sumitomo 
F-70H helium compressor fE3l1 . An outer vacuum chamber at a pressure of 10“^ mbar isolates 
the cryogenic refrigerator from ambient room temperature. The purification of the argon gas is 
achieved with two purification cartridges OXISORB-W mounted in parallel, which reduce the O 2 
and H 2 O levels to < 5 and < 30 ppb respectively fE^ . The temperature is kept between the so¬ 
lidification and the boiling points of argon (85 K and 87 K respectively at 1 atm) by four resistive 
heaters located on the first stage cooling station of the cryocooler. Platinum sensors are connected 
along the coil pipe to monitor the cooling power. A PID controller driven by a Lab View program 
reads out and regulates the temperature of the cold head at the desired value with a precision of 
about 10 mK. The current provided to the heaters is ensured by a TTi TSX 35 10 programmable 
PSU based on the GPID interface fE5l1 . The typical power produced by the heaters to maintain the 
temperature for liquefaction is about 50 W. 

The neutrons scattered off the LAr cell are detected by a 5" organic liquid scintillator counter 
(LSC) manufactured by SCIONIX fEHl . The scattering angle © is varied by rotating the LSC 
around the LAr cell on a 1 m long arm. The LSC is shielded against direct neutrons from the 
source by a 10 cm thick polyethylene absorber. Events are recorded by a coincidence measurement 
of the LAr and LSC signals. To obtain the ratio Leff, equation (pn|), the light output produced 
by the argon nucleus emitted under the corresponding recoil angle is compared to the reference 
electron recoil light yield. 



Fast Amp LEO discri Logic unit Timer 


Figure 4. Diagram of the data acquisition system. 


Figure ^ shows the schematic of the data acquisition system. The analog signals from the 
two LAr PMTs are split (not shown in the figure) and fed into the low and high gain inputs of a 
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10-bit digitizer (LeCroy Oscilloscope WavePro 735Zi DSO). The low and the high gain signals 
are matched later offline. This method increases the dynamic range so that signals from single 
photoelectrons and from a-emitters are collected with high sensitivity without ADC saturation 
lETII . A cross check is performed by comparing the pulse height obtained by this method with the 
pulse height from several y-sources. The accuracy of the signal matching between the high and 
low gain is ^2%. The data are collected with 5 000 points at IGS/s with 500 ns pre-sample for 
pedestal subtraction. Data are stored on a hard disk and then analyzed offline. The signals from 
each PMT are fed into a CAEN N979 xlO amplifier before being fed into leading edge CAEN 
N840 discriminators where the thresholds have been set to 0.5 photoelectrons. The logic signals 
are then connected to a CAEN N405 unit for the coincidence between the two PMTs. 

Neutron signals from the ESC are processed by an analog pulse shape discriminator module 
Mesytec MPD4 lESlI . This module determines the fraction of fast (<20 ns) component of the 
scintillation light. Eor electrons the prompt fraction Fp (prompt over total integrated pulse height) 
is about 0.95, while for proton recoils Fp < 0.8 (figure g). 



Integrated pulse height [nVs] 

Figure 5. Prompt fraction Fp of scintillation light in the ESC versus integrated pulse height. The upper 
band is produced by electronic recoils from bremsstrahlung (generated by the neutron generator) and cosmic 
muons, the lower band by proton recoils from neutron scattering. 


The energy calibration of the ESC is performed with radioactive y-sources. The energy resolu¬ 
tion is determined from the signals induced by backscattered photons from ^^^Cs or ^^Na sources, 
which are detected by a small BaE crystal lE9l1 . The time-of-flight between the EAr cell and the 
ESC is measured to select elastic events. The ESC signal is fed into a linear fan in/out module 
which feeds the timer and logic units. The time-of-flight between the EAr cell and the ESC is 
calculated from the time difference of the TAC output of the MPD4 module and the signal from the 
EAr cell. The time-of-flight is calibrated by the coincidence of the 511 keV back-to-back y-rays 
from a ^^Na source located at mid-distance between the EAr cell and the ESC. Events are accepted 
within a time delay of 200 ns between the EAr and ESC signals. This trigger setting was used 
during part of the data taking (TR 2 data sample, see table |T| below). Eor the TRi data sample taken 
earlier a programmable trigger logic in the oscilloscope was used for the coincidence between the 
EAr and ESC signals. 
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3. Calibrations with y-sources 


3.1 Pulse height analysis 

The light yield from the LAr cell is the ratio of collected to produced light intensities, which 
strongly depends on factors such as geometry, photomultiplier efficiency, wavelength shifter, etc. 
The energy calibration in LAr is performed by using the 60 keV photo-absorption peak from an 
external Am source. The source is placed at a distance of 10 cm from the vessel to reduce pile- 
up events. The light yield has to be corrected for the finite integration time. As mentioned earlier, 
nuclear recoils contribute mostly to the fast component (Ti), while the contribution from electrons 
to the slow component (T 2 ) is larger. Therefore the light yield has to be corrected first for losses 
due to impurities in LAr which are reducing the lifetime of the slow component fEOl . 



Figure 6. Left: logarithmic binning method showing the contributions of the populations A, B and C. Right: 
pulse height distribution from the 60 keV 7 -source (in photoelectrons, p.e.) for the raw data (black curve) 
and after correction with the log binning method (red curve). 


We have developed a new technique to correct for the reduced light yield caused by argon 
impurities. This event-by-event method consists in dividing the pulse shape in three time regions 1 \, 
I 2 , and L, as sketched in figure ^ (left). Time / = 0 corresponds to -20 ns below the maximum pulse 
height. The first interval (A) corresponds to 4 ti. The integral over A contains the contributions 
A from the fast component and B from the slow component. A third contribution C is present in 
the interval I 2 , between 40 ns and 100 ns. This component has been observed earlier in argon 
and is reported to stem from the wavelength shifter fCTl . The rest of the pulse shape is divided 
logarithmically into n = 5 time bins Li from /q to defined as 

U — t{)-\-'l2’\n( ——-: I (3.1) 

with / = 1...5 and where T 2 is the measured decay time of the slow component. Thus the expected 
number of events in each bin L/ is constant. Impurities affect the late regions of the light pulse, 
hence affect only the slope of the slow component B but not its amplitude A/q at / = 0. The measured 
number of entries in Li can be obtained by building the arithmetic mean of the five bins. The con¬ 
tribution Bcorr from the slow component, corrected for impurities, is then obtained by integrating 
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N{t) with the time constant T 2 = 1.6 /is which corresponds to pure LAr. Once Bcorr is obtained, C is 
found from the h integral and, finally, A from the h integral after subtracting the contribution from 
Bcorr and C. Details can be found in flTlII . Let us define the fraction of prompt light determined 
with the logarithmic binning by the ratio 

A + Bcorr+ C' 

Figure ^ (right) shows the raw energy distribution from the Am source for runs taken under 
various purity levels. Mixing data with different purities leads to a spread of the deposited energy 
distribution (black curve). However, applying the logarithmic binning method to correct the spec¬ 
trum leads to a decrease of the width of the distribution and shifts the distribution towards higher 
photoelectron numbers. 

Light yield calibrations were performed periodically during data taking. The corrected aver¬ 
age light yields are 3.75 ±0.08 p.e./keVg^ and 3.39 ±0.07 p.e./keVg^ for the TRi and TR 2 data 
respectively flTHl . The slightly different light yields are due to the different thicknesses of wave¬ 
length shifter covering the surfaces of the PMTs. The long term stability of the light yield is shown 
in figure 0 for the TRi data. Further external radioactive sources were employed to cover a wide 
energy range: 122 keV photoelectron from ^^Co, 341 keV, 478 keV and 1060 keV Compton edges 
from ^^Na, ^^^Cs and ^^Na respectively. Figure ^ illustrates the perfect linearity of the system. 
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Figure 7. Light yield as a function of time for the TRi data sample. 



Figure 8. Light yield from LAr measured with various 7 -sources (in arbitrary units). In the present work 
the calibration is obtained from the light yield of the 60 keV Am photopeak. 
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The response to very low energy electrons was also measured with a source (half-life 
of 86 days) of about 180 kBq which was introduced into the system. A ^^Rb trap was directly 
connected to the gas system. The daughter metastable ^^^Kr nuclide (half-life of 1.8 h) illuminates 
the central part of the LAr cell and diffuses photons uniformly. The measured light yield at 32.1 
keV is obtained as 3.68 ±0.15 p.e./keV^^ in good agreement with results from the other sources. 

3.2 Trigger efficiency 

Understanding the low energy roll-off of the trigger efficiency is essential to determine Le/f, espe¬ 
cially for low energy recoils. We have studied the trigger efficiency with the 511 keV annihilation 
y-rays from a ^^Na positron source. The energy distribution of the Compton scattered electron is 
linear at low energy, as illustrated by the red dashed box in figure ^ (left). As a comparison figure 
0 (right) shows the measured electron recoil spectrum in LAr. 



Figure 9. Left: theoretical electron recoil spectrum from the 511 keV annihilation y of a ^^Na positron 
source (computed with the Klein-Nishina formula). The red dashed box shows the linear part of the spectrum 
at low energy. Right: measured integrated pulse height v of the prompt light in photoelectrons (p.e.) fitted 
with the function ([3.3|) (red curve). The inset shows the low energy region. 

The low energy part is fitted with the following model function: 

^(x) =5-[(D±^-x)±L(a,Ai,cj)]. (3.3) 

The function 

5=(l-e-^)" (3.4) 

describes the trigger efficiency, where x is the pulse height in photoelectrons. The term {D + E -x) 
describes the linear shape of the spectrum at low recoil energies. The parameters D and E are 
determined by the fit. The Landau distribution takes the Cerenkov light from the y- 

ray interaction in the PMT glass into account, where a is a normalization parameter, the most 
probable value and a a scale parameter (adapted from the CERNLIB routine GllO denlan). Most 
of the Cerenkov events contribute to the prompt light. 

The fit is performed in two steps. First, the linear part of the spectrum is fitted. Second, the 
parameters of the linear function and the Landau function are fixed and the parameters b and c are 
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Figure 10. Trigger efficiency S for the source measurement as a function of integrated pulse height x 
for the two trigger settings. 

deternfined. The measured trigger efficiency S is displayed in figure [T^ by the green curve for the 
S 2 data. The efficiency is ^100% at 8 photoelectrons (p.e.) , ^96% at 4 p.e., ^50% at 2 p.e. and 
^5% at 1 p.e. The roll-off of the trigger efficiency for the TRi data that used the progranunable 
trigger logic in the oscilloscope is shown by the blue curve. The efficiency is lower, reaching ^95% 
at 8 p.e., ^55% at 4 p.e. and ^6% at 2 p.e. 

4. Neutron data 

The neutron generator was operated most of the time at 80 kV to keep the accidental background 
from bremsstrahlung generated by the neutron generator at an acceptable level. The corresponding 
neutron flux is 6 x 10^ neutrons/s into 471. The scattered rate from the sensitive LAr volume, 
recorded by the LSC, is 1 neutron/min. The background rate is estimated to be 5/min, caused 
mainly by diffusely scattered neutrons and X-rays from bremsstrahlung (^3-4/nfin), and to a lesser 
extent by cosmic muons saturating the MPD4 analog pulse shape discriminator. 
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Trigger 
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25 

11.5 

TRi 

269 

TR 2 

382 

30 

16.4 

TRi 

284 

TR 2 
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40 

28.5 

TRi 

379 

TR 2 

219 

50 

43.4 

TRi 

291 

- 

- 

60 

60.5 

TRi 

288 

- 

- 

90 

119.5 

TRi 

196 

- 
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Table 1. Angular settings and corresponding argon recoil energies. The trigger type (see text) and running 
time are also given. 


Neutron scattering data are measured at six angles corresponding to argon recoil energies 
between 11 and 120 keV (table |T|). The data are taken with the two different triggers TRi and TR 2 . 
The raw data are first corrected for impurities and for the trigger efficiency roll-off described in 
section E3- 
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The discriminator vlb (23), which is based on the logarithmic binning method (and thus differs 
from the prompt fraction Fp), suppresses background from accidental coincidences with X-rays, 7 - 
rays and cosmic muons. As an example, figure shows a scatterplot of vlb in LAr as a function of 
integrated pulse height for the © = 25^^ data sample. Events above the green line are mostly neutron 
induced, those below the line stem mostly from background. 



Figure 11. Scatterplot of the riB discriminator versus integrated pulse height for 0 = 25®. The green line 
shows the cut applied to select elastic neutron scattering (red region) and reduce background. 


Figure (left) shows a scatterplot of the prompt light fraction Fp (defined in section ^ versus 
time-of-flight for © = 25^^. The event distributions for the data taken at the larger scattering angles 
are similar. The time-of-flight projection is shown in figure (right). Events in red are nuclear 
recoils, while events in blue are mostly due to background (lying below the green line in figure [TTI) . 
Since the LAr and the ESC triggers are not correlated the accidental background is uniform and 
appears as a flat component in the time-of-flight spectrum. 



Time of flight [ns] 

Figure 12. Prompt fraction versus time-of-flight (left) and corresponding projection on the time-of-flight 
axis (right) for 0 = 25®. The events in red correspond to nuclear recoils selected by the vlb cut, while the 
anticut selects the electron-recoils (events in blue). The events in green correspond to nuclear recoils at 
energies above 120 photoelectrons. The elastic neutron scattering events selected to calculate L^/f lie in the 
pink window. The vertical dashed lines delimit the intervals used to evaluate the contribution from accidental 
events. 
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The most prominent red peak is due to elastically scattered neutrons at the expected time- 
of-flight of about 43 ns. These events are also associated with large values of Fp and therefore 
contribute mostly to the prompt light. An accumulation of events is observed around 80 ns, corre¬ 
sponding to inelastically scattered neutrons in LAr. The broad tail of this distribution is mainly due 
to neutron scattering outside the active volume and inelastically scattered neutrons. The two peaks 
around zero time-of-flight are due to cosmic muons flying between the LAr cell and the LSC, the 
peak at -3.1 ns corresponding to muon from the LSC to the LAr cell and the one at +3.1 ns to those 
flying in the opposite direction. In the electronic (blue) recoil band an accumulation of events is 
observed above the elastic peak. They are correlated in time with neutrons scattered inelastically 
in the target but outside the sensitive volume, while the ensuing emitted photons are converted in 
the active volume. The events shown in dark green are nuclear recoil events with energy deposits 
above 120 photoelectrons, mostly due to multiple neutron scattering. The accumulation of events 
at -50 ns is due to an electronics artifact at the edge of the trigger time window. 

Two cuts are applied to select elastic neutron scattering. A time-of-flight cut window is applied 
between -2 ns and +6 ns of the most probable value for the elastic peak (pink area in figure [T2| 
corresponding to 41 - 48 ns at 0 = 25^). The cut by the ris discriminator (green line in figure 
OD removes efficiently accidental coincidences with X-rays, y-rays, inelastic scatters and cosmic 
muons. Figure shows the contribution to the fast light for the selected data. The recoil energy 
distribution shown in this figure will be fitted to the Monte Carlo data, as described in section 



Integrated pulse height [p.e.] 

Figure 13. Nuclear recoil energy distributions at 0 = 25® (in photoelectrons). The plot shows the contri¬ 
bution to the fast light. The accidental spectrum is shown in blue. The vertical red dashed line indicates the 
pulse height corresponding to 90% trigger efficiency. 


To estimate the contribution from residual accidentals we choose events lying far from the 
elastic peak in the two intervals delimited by the vertical dashed lines in figure [^. The background 
spectra below and above the elastic peak are found to be compatible. The contribution from acci¬ 
dentals is negligible, as shown by the blue histogram in figure |^. However, neutrons that scatter 
off various materials outside the LAr cell before or after interacting in the LAr sensitive volume 
also contribute background. This external background cannot be removed from the data and is 
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therefore taken into account by Monte-Carlo simulation. 


5. Monte-Carlo simulation 

The detector response was simulated extensively with GEANT4. The Monte-Carlo simulation 
takes into account the neutron generator assembly with its shielding and collimator, the cryostat 
and the various components of the LAr cell such as the vacuum vessel, the TPB reflector foil, the 
photomultipliers, the support mechanics, and the LSC. A photograph of the LAr cell and a drawing 
of the simulated one are compared in figure 0 . 


CH2 holder 
PMT 

Al cylinder 
LAr volume 

PMT 



Figure 14. Photograph of the LAr cell and PMT (left) and drawing of its simulated counterpart (right). 


In the simulation monoenergetic neutrons of 2.45 MeV are emitted isotropically into 47r. A 
sample of 100 million interactions in the active LAr volume was generated at each measured angle. 
Figure (left) shows for 0 = 25^ the various simulated contributions to the nuclear recoil energy 
spectrum (left) and to the time-of-flight distribution (right) prior to smearing by the experimental 
resolution. A clear elastic peak is observed at the expected recoil energy, as well as the exponen¬ 
tially decreasing material background tail (black histogram). Neutrons that scatter elastically only 
once in the active LAr volume are shown by the blue histogram. The dashed red histogram corre¬ 
sponds to single scattering in the active LAr volume, preceded or followed by scattering elsewhere 
in the apparatus (external background). The multiple scattering events in LAr are shown in green. 
Inelastically scattered neutrons in LAr are represented by the dashed orange histogram. The energy 
deposits are shifted to higher values by the conversion of the 7 -rays emitted during the de-excitation 
of the argon atoms. 

Figure (right) shows the corresponding time-of-flight distributions at © = 25^^. The position 
in time of the elastic peak is consistent with data. A comparison with the data of figure shows 
that nearly all inelastic scatters are removed by the time-of-flight cut of 41- 48 ns. Their time-of- 
flight distribution (enhancement around 70 ns) is also in good agreement with data. 

The position of the single elastic peak and its resolution depend on the sizes of the neutron 
source and on the solid angle acceptances of the LSC and LAr cell. Uncertainties on recoil energies 
are determined by fitting a Gaussian to the single elastic peak. Before applying the time-of-flight 
cut to the simulated data the finite time resolution is taken into account by convoluting a Gaussian 
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Figure 15. Simulated recoil energy distribution at 0 = 25® (left) and corresponding time-of-flight spec¬ 
trum (right) from the contributions listed in the inset: single elastic scattering (with and without external 
scattering), multiple scattering and inelastic scattering in the LAr sensitive volume. 


with the experimental rms resolution of 2.88 ns. The recoil energy distribution at © = 25® is shown 
in figure after the time-of-flight cut. The black histogram will be used to fit to the data of figure 

o. 



Figure 16. Simulated recoil energy distribution at 0 = 25® after the time-of-flight cut. The color codes are 
the same as in figure |^. 


The various contributions to the nuclear recoil spectrum after the time-of-flight cut are listed 
in table 0 for the six scattering angles. Below 40® single elastic scattering off argon is the dominant 
contribution to the recoil energy spectra. The contribution from multiple scattering in LAr is small 
(^10%) thanks to the small dimension of the LAr cell compared to the elastic neutron mean free 
path. The decreasing contribution of elastic events with increasing scattering angle is due to the 
angular dependence of the differential elastic cross section for neutron-argon interaction at 2.45 
MeV which reaches a minimum around 90° lE2]. 
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© 

Single 

Single elastic 

Multiple 

Inelastic 


elastic 

+ external 

elastic 


[“] 

[%] 

[%] 

[%] 

[%] 

25 

67.7 

23.5 

8.8 

0.02 

30 

65.3 

25.1 

9.6 

0.02 

40 

54.8 

31.4 

13.7 

0.06 

50 

49.7 

34.3 

16.0 

0.04 

60 

40.9 

31.8 

20.5 

0.04 

90 

29.1 

45.1 

25.6 

0.18 


Table 2. Rms contributions to the recoil spectrum after time-of-flight cut. 


6. Data analysis 


6.1 ;t^-fits 


The method developed in this work to extract Leff follows the one described in [HBL13] for liquid 
xenon. We have performed a ;^^-niinimization to determine the relative light efficiency Leff ( ^T| ) 
as a function of recoil energy Tnr at six scattering angles. 

To obtain the simulated spectrum in photoelectrons the energies in keV are first multiplied by 
the free parameter Leff under tost and then by the light yield Ly measured with the Am source. 
The number of photoelectrons Np^e. is allowed to fluctuate according to a Poisson distribution and 
the energy resolution is described by the rms deviation of a Gauss distribution defined as R-s^Np^e., 
where R is the second free parameter under x'^ test. The gain fluctuation of the PMTs is taken 
into account by convoluting the nuclear recoil energy spectrum in photoelectrons with a Gaussian 
distribution with rms deviation Rp e. -\/Np^e.- The parameter Rp e. is determined from the resolution 
of the measured PMT single photoelectron distribution and found to be Rp ^. = 0.4 llTHl . The trigger 
efficiency (figure [T^) is then taken into account to obtain the simulated recoil energy spectrum. For 
each scattering angle the x^ is computed as follows: 


XHLeff,R) 


N 


L 


[hi-hMc,i{Leff,R)]^ 


( 6 . 1 ) 


where Leff = Leff(Tnr) and the resolution R = R(Tnr) are the free parameters. The events are 
divided into N bins expressed in photoelectrons. The number of events in the measured and the 
simulated bins i are labelled hi and Hmcj respectively. The corresponding statistical uncertainties 
are denoted by Ot and cJmc,/- The fit range is chosen separately for each angle in such a way as to 
avoid the bias induced at low energy by the trigger inefficiency (which underestimates Leff) and 
by the high energy tail (which overestimates Leff). Fits are performed for different photoelectron 
ranges and the suitable fit range is selected as the one associated with stable values of Leff. For 
example, at 25^ the range fitted was from 2 to 38 photoelectrons (see figure The x^ surfaces 
are illustrated in figure (left) as a function of R and Leff for © = 25^ (S 2 data), 40^ (S 2 data) and 
90^ (TRi data). 

Following the method described in llT9l1 a rotated paraboloid is fitted to the x^ generated 
around the minimum value for each of the six scattering angles &j. This procedure to extract the 
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Figure 17. Left: maps as a function of and R for three scattering angles. The location of the 

minimum is shown by a dot. The curves are the la and 2a contour lines. Right: corresponding Monte- 
Carlo fits (red histograms) to the energy recoil spectra (black data points) with rms errors (grey zones). The 
dotted blue lines delimit the fitted ranges. 


final values for L^f/j and Rj (j = 1...6) averages out small variations on the x^ surface. The rotated 
X'^ is defined as 


X\Leff,R)=xlin,j + 
+ 


i^eff ^effj) COS (Oj {R Rj) sincOj 


i^eff ^e//,;') sinrOj {R Rj) cos CO J 
Bi 


( 6 . 2 ) 
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where the paraboloid rotation angle (Oj and the parameters Aj and Bj are determined by the fit, 
together with the minimum Zminj- "*"he location of Xminj indicated by the black dots in figure 
(left). The la and 2a contours, determined by the values Xmin + 2-3 and Xmin + respectively, 
are also shown in figure [T^ (left). 

Figure [T^ (right) shows the measured nuclear recoil distributions (black dots) and the Monte- 
Carlo fits (red histograms). The grey zones show the statistical uncertainties on the fit, determined 
by the la contours in figure |T7| (left). The blue dashed lines correspond to the fit ranges used to 
compute the x'^- The results of the fits are given in table || with statistical errors only. 


© 

["] 

^eff 

^eff 

(stat) 

R 

A/? 


d.o.f 

25 

0.386 

0.0185 

2.26 

0.29 

18.0 

16 

30 

0.305 

0.0142 

1.84 

0.28 

19.8 

18 

40 

0.285 

0.0127 

1.86 

0.30 

20.9 

19 

50 

0.294 

0.0159 

2.79 

0.42 

36.8 

24 

60 

0.283 

0.0167 

3.43 

0.72 

33.6 

20 

90 

0.301 

0.0177 

3.49 

0.74 

15.4 

10 


Table 3. Results of the fits. Listed are the optimum values of the scintillation efficiency L^ff and resolution 
R as a function of scattering angle 0. The rms errors are statistical only. The x^ values and numbers of 
degrees of freedom (d.o.f) refer to the least square fit (|0|) . 


6.2 Systematic errors 

The reliability of the procedure to extract Le/f by including the additional fit parameter R has been 
studied in flT^II for liquid xenon. The analysis was cross-checked by describing the energy depen¬ 
dence of the resolution with model functions without free parameters. No significant deviation 
was observed between the values of Leff obtained with this method and those determined with the 
resolution left as a free parameter. A similar test is performed in our work: the measured nuclear 
recoil spectra with the elastic peak above trigger efficiency roll-off are fitted with Gaussian func¬ 
tions to obtain the experimental resolution. The simulated contribution of the energy spread due 
to the finite size of the detector is then subtracted quadratically from the measured resolutions. No 
significant deviation is observed on Leff between the two methods (figure [rS| ). 

The following sources of systematic uncertainties are taken into account. The uncertainty in 
the light yield ALy measured with the ^"^^Am source is determined by the error on the calibration 
constant, estimated to ±2.1% for the TRi and ±1.6% for the TR 2 data. The uncertainty AT rig on 
the trigger efficiency is determined by changing the values of the parameters b and c in The 
error A© on the angular setting of the LSC is estimated to be ± 0.5^. The systematic uncertainties 
are listed in table 0 together with our final results. The main contributions to the systematic error 
on Leff stem from the calibration of the americium source and from the trigger efficiency roll-off. 
The latter large value at 50^ is due to the trigger settings for the TRi data, for which the effect of 
the roll-off is more pronounced (see figure |^. To estimate the systematic error introduced by the 
uncertainty on the lifetime of the slow component we have varied T 2 between 1.5 and 1.6 /is. This 
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Figure 18. Plot showing how Leff varies when using fixed resolutions R (blue circles) instead of free fit 
parameters (black rectangles). 


leads to an additional (negligible) systematic error of typically 0.003 on Leff. The uncertainties in 
nuclear recoil energies are determined by the angular acceptance of the LSC and are estimated by 
Monte-Carlo simulation. The total uncertainties on Leff are obtained by summing the statistical 
errors (table ||) and systematic errors quadratically. 


©n 

T„r[keV] 

^eff 

^eff 

ALy 

AT rig 

A© 

25 

11.5±2.8 

0.386 

0.032 

0.0068 

0.0252 

0.0028 

30 

16.4±3.9 

0.305 

0.020 

0.0054 

0.0126 

0.0010 

40 

28.5±5.2 

0.285 

0.014 

0.0051 

0.0013 

0.0007 

50 

43.4±6.0 

0.294 

0.023 

0.0064 

0.0154 

0.0006 

60 

60.5±7.1 

0.283 

0.018 

0.0062 

0.0022 

0.0004 

90 

119.5±7.7 

0.301 

0.019 

0.0064 

0.0000 

0.0000 


Table 4. Scintillation efficiency Leff and rms error AL^// as a function of scattering angle. The last three 
columns list the rms systematic errors to ^eff- 


6.3 Discussion 

The final results for the relative scintillation efficiency from table ^ are plotted in figure together 
with previous measurements. Above :^25 keV recoil energy Leff is constant with a mean value 
{Leff) = 0.30 zb 0.02. This result is in fair agreement, although somewhat higher, than that from 
MicroCLEAN which reports an average value of 0.25 zb 0.02 above 20 keV flT^ . Our simulations 
are in excellent agreement with data and we stress here the importance of modelling accurately in¬ 
teractions outside the active volume, in particular for large active volumes (3.14 ^ in firai compared 
to 0.2 I in the present work). 

Below 20 keV our data show the upturn in Leff already reported by MicroCLEAN flTHl but 
not observed by SCENE flTTlI . Leff is usually predicted to decrease at low energy due to the com¬ 
bination of energy loss (described by Lindhard’s theory) and scintillation quenching (described by 
Birk’s law) [BUI . This prediction is shown by the dotted curve in figure |^. However, based on the 
work described in [B5l1 the NEST group has recently applied on argon their model for liquid xenon 
[BHi . The authors developed a parametrization to compare with liquid xenon data that reproduce 
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Figure 19. Relative scintillation L^ff in liquid argon as a function of nuclear recoil energy from the present 
experiment (red squares), compared to the previous measurements by MicroCLEAN (blue triangles) flT^ . 
SCENE (black diamonds) at zero electric field [O] and WARP (open circle) O. The dotted and dashed 
curves show the predictions from the theoretical models by Mei and NEST . 


the upturn observed here and by MicroCLEAN. The upturn could for instance be due to an increas¬ 
ing exciton-ion ratio at low energies, which would lead to an increasing light yield since it takes 
less energy to excite than to ionize fETlI . The NEST prediction extrapolated to LAr fCTi is shown 
by the dashed curve in figure |^. 

7. Conclusions 

Summarizing, the knowledge of the scintillation efficiency Leff at low nuclear recoil energies is 
important for direct dark matter searches using noble liquids. We have measured Leff for nu¬ 
clear recoils, relative to electrons, between 11 keV and 120 keV in liquid argon. Single elastic 
neutron-argon scattering dominates thanks to the small active volume of our argon cell and the 
well collimated neutron beam, while the contamination from neutron inelastic scattering is negligi¬ 
ble. A minimization is performed leaving Leff and the energy resolution as free parameters. The 
extraction of Leff is challenging at very low energy where systematic uncertainties on Leff are de¬ 
termined by the inefficiency of the trigger. At higher energy the uncertainties on Leff are dominated 
by statistical errors. The scintillation efficiency is constant with mean value {Leff) = 0.30 ±0.020 
between 16 keV and 120 keV. The results below 20 keV confirm the energy upturn reported earlier. 
The increasing value of Leff will enhance the detection efficiency for low mass WIMPs and be 
beneficial to dark matter searches using LAr. 
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